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Abstract. Libraries of stellar spectra are fundamental tools for the study of stellar populations and both empirical 
and synthetic libraries have been used for this purpose. In this paper, a new library of high resolution synthetic 
spectra is presented, ranging from the near-ultraviolet (300nm) to the near-infrared (1.8/jm). The library spans 
all the stellar types that are relevant to the integrated light of old and intermediate-age stellar populations in 
the involved spectral region (spectral types F through M and all luminosity classes). The grid was computed for 
metallicities ranging from [Fe/H] = -2.5 to +0.5, including both solar and a-enhanced ([a/Fe] = 0.4) chemical 
compositions. The synthetic spectra are a good match to observations of stars throughout the stellar parameter 
space encompassed by the library and over the whole spectral region covered by the computations. 
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1. Introduction 

Libraries of stellar spectra are one of the main ingredi- 
ents of stellar population synthesis models (e.g. Bruzual 
& Chariot 1993, 2003; Cervino & Mas-Hesse 1994; Fioc 
& Rocca-Volmerange 1997; Leitherer et al. 1999; Vazdekis 
1999; Buzzoni 2002; Schulz et al. 2002; Jimenez et al. 2004; 
Gonzalez Delgado et al. 2005; Maraston 2005; Schiavon 
2005) and both empirical and theoretical libraries have 
improved dramatically in recent years. The first empir- 
ical libraries (e.g. Gunn & Stryker 1983; Jacoby et al. 
1984) used in stellar population synthesis work were re- 
stricted to a relatively small number of stars with uncer- 
tain atmospheric parameters. The quality of empirical li- 
braries has been refined along the years (e.g. Pickles 1998; 
Jones 1999) and recently, a major improvement has been 
achieved by Le Borgne et al. (2003; STELIB) and Valdes 
et al. (2004; Indo-US), which provide high S/N, medium 
resolution (down to FWHM ~ lA) and a good coverage 
of the color-magnitude diagram. 

Amongst the synthetic libraries, perhaps the most 
widely used is the flux distribution predicted from the 
Kurucz (1993) model atmospheres. Lejeune et al. (1997, 
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1998) and Westera et al. (2002) extended this library to 
include spectra of M stars computed from model atmo- 
spheres from Fluks et al. (1994), Bessell et al. (1989), 
Bessell (1991) and Allard & Hauschildt (1995). Moreover, 
they calibrated the flux distribution of the synthetic spec- 
tra in order to match the colors of observed stars (BaSeL 
library). However, the spectral resolution of the BaSeL li- 
brary is limited to the sampling of the model atmosphere 
grid (~ 20 A) which is by far lower than the modern ob- 
served spectra of both individual stars and integrated stel- 
lar populations. More recently, resolution ceased to be a 
limitation, with the publication of high resolution spec- 
tral libraries by Chavez et al. (1997), Murphy & Meiksin 

(2004) , Martins et al. (2005) and Rodrigues-Merino et al 

(2005) . However, these theoretical libraries still have a 
more limited wavelength coverage than that of the pre- 
vious low-resolution synthetic libraries. 

The choice of using either an empirical or a synthetic 
library in stellar population models is a subject of debate 
in the literature. One disadvantage of synthetic libraries 
is that they rely on model atmospheres, which are sub- 
ject to systematic uncertainties. Besides, computing a re- 
liable high-resolution synthetic spectral library for a large 
range of stellar parameters and in a wide spectral region 
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is a very challenging task, because it requires building an 
extensive and reliable list of atomic and molecular line 
opacities which are needed for an accurate reproduction 
of high- resolution spectra of real stars. On the other hand, 
synthetic spectral libraries overcome limitations of empir- 
ical libraries, the most important being their inability to 
extrapolate to abundance patterns that differ from that 
of the library stars, which are from the solar neighbour- 
hood or, in some cases, that of the Magellanic Clouds. 
Therefore, models based on empirical libraries cannot re- 
produce the integrated spectra of systems that have un- 
dergone star formation histories that are different from 
that of local systems. 

The first compelling evidence that models based on 
empirical libraries have difficulty reproducing the inte- 
grated light of extra-galactic populations was presented 
by Worthey et al. (1992). These authors showed that sin- 
gle stellar population synthesis models for the Lick/IDS 
indices (Burstein et al. 1984; Worthey et al. 1994; Trager 
et al. 1998) cannot reproduce the indices measured in gi- 
ant elliptical galaxies, thus indicating that these systems 
are overabundant in a-elements relative to the Sun. This 
happens because, by construction, the abundance pattern 
of models based on empirical libraries is dictated by that 
of the library stars, which mirrors the abundance pattern 
of the solar neighborhood (but sec discussion in Schiavon 
2005). Therefore, such models are bound to follow the 
relation between [Fe/H] and [a/Fe] that is characteristic 
of the solar neighborhood where [a/Fe] is above solar at 
low metallicities and around solar in the high-mctallicity 
regime (e.g. McWilliam 1997). As a result, the models can- 
not match the integrated light of metal-rich, a-enhanced 
systems, such as giant ellipticals. 

This deficiency could be in principle partially cured 
empirically, by including spectra of bulge stars in the spec- 
tral libraries. The population of the bulge of the Galaxy 
has an a-enhanced abundance pattern even at high [Fc/H] 
(e.g. McWilliam & Rich 1994, 2004; Zoccali et al. 2004). 
However, their inclusion in sufficient numbers is hampered 
by their faintness and severe reddening, so that very few 
(or none) of these stars are included in empirical libraries. 

For the case of the Lick/IDS indices, models account- 
ing for variable a-enhancement have been developed (e.g. 
Trager et al. 2000; Thomas et al. 2003a; Proctor ct al. 
2004; Tantalo et al. 2004; Mendes de Oliveira et al. 2005; 
Schiavon 2005). Even though the Lick/IDS indices have 
already proven their critical importance for the under- 
standing of the stellar populations, it is crucial to develop 
models that can explore the huge amount of information 
contained in the whole spectrum, since methods fitting 
the entire spectrum are now feasible (Heavens et al. 2004; 
Pantcr ct al. 2004). 

Currently, the only way of producing models with 
abundance patterns that differ from that of the stars in the 
solar neighbourhood is through the adoption of synthetic 
stellar spectra. Although grids of high-resolution synthetic 
spectra have recently been computed, a wide grid with an 
a-enhanced mixture is still lacking. To our knowledge, so 



far two grids of synthetic spectra with a-enhanced mix- 
tures were available: Barbuy et al. (2003) published a grid 
in the wavelength range 4600 - 5600 A aiming at the study 
of the effect of a-enhancement on the Lick/IDS indices 
Mg2, Fe5270 and Fe5335; and Zwitter et al. (2004) pub- 
lished a grid covering 7653 - 8747 A, which is aimed to be 
used with the spectroscopic surveys by RAVE and GAIA 
(Steinmetz 2003; Katz 2004). 

In this work we present for the first time a grid of 
synthetic spectra covering solar and a-enhanced mixtures 
in a wide baseline. The grid is aimed at applications in 
stellar population synthesis, specifically to the study of 
old and intermediate-aged stellar populations. Thus, the 
grid covers the effective temperatures and surface gravities 
suitable for these systems. 

We would like to emphasize, however, that applications 
of this grid are not limited to stellar population studies. 
The spectral library should also be a very useful tool for 
atmospheric parameter determination work, in particular 
as input in automatic procedures of deriving atmospheric 
parameters of a large number of observed stars (e.g. Cayrel 
et al. 1991a; Franchini et al. 2004; Willemsen et al. 2005; 
Girard & Soubiran 2005; Valenti & Fischer 2005). 

This paper is organized as follows. In Section 2 we 
briefly describe the computation of the synthetic spectra 
and the production of the atomic and molecular line lists. 
In Section 3 the library is presented and aspects of the rel- 
ative contribution of the chemical species are investigated. 
A Summary is given in Section 4. 

2. Calculations 

The spectra were computed using the code PFANT. The 
first version of this code named FANTOM was developed 
by Spite (1967) for the calculation of atomic lines. Barbuy 
(1982) included the calculation of molecular lines, imple- 
menting the dissociative equilibrium by Tsuji (1973) and 
molecular line computations as described in Cayrel et al. 
(1991b). It has been further improved for calculations of 
large wavelength coverage and inclusion of Hydrogen lines 
as described in Barbuy et al. (2003). Ten hydrogen lines 
from H w to H a are considered, through a revised version 
of the code presented in Praderie (1967). Higher order 
Balmer lines and Paschen and Brackett series are not yet 
considered in the present version of the code. Given a stel- 
lar model atmosphere and lists of atomic and molecular 
lines, the code computes a synthetic spectrum assuming 
local thermodynamic equilibrium. 

The model atmospheres adopted were those presented 
by Castelli & Kurucz (2003), which account for a- 
enhanced compositions. The solar abundances adopted 
are those from Grevesse & Sauval (1998), which are con- 
sistent with the abundances used in the model atmo- 
spheres. These models differ from the previous ATLAS9 
grid (Kurucz 1993) mainly because they include new opac- 
ity distribution functions, new molecular opacities and do 
not include overshooting. 
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Fig. 1. Upper panel: comparison of the observed solar spectrum (thin line) and synthetic spectrum (thick line). Lower 
panel: comparison of Arcturus observed spectrum (thin line) and synthetic spectrum (thick line) . Adopted parameters 
are given in the text. 



These models have mixing length to pressure scale 
height ratio of 2.0. A lower value (I /Hp — 0.5) was 
suggested by Fuhrmann et al. (1993) and Van 't Veer- 
Menneret & Megessier (1996) as being more suitable to 
reproduce the profile of Balmer lines (as was adopted in 
the models in Barbuy et al. 2003). But it is important to 
note that the Balmer lines are not well reproduced also 
because their cores are affected by the chromosphere and 
by departures from local thermodynamic equilibrium. 

2.1. Atomic line opacities 

Accurate line opacities are in the very heart of a successful 
synthetic spectrum computation. The line list employed is 
an updated version of the one presented in Barbuy et al. 
(2003), where oscillator strengths (log gfs) and damping 
constants were revised. 

The line list was extended up to 1.8/im with the in- 
clusion of the data from Melendez & Barbuy (1999). The 
most recent oscillator strengths listed in the latest elec- 



tronic version of the NIST database (Reader et al. 2002) 
were adopted. These values account for 10 to 50% of the 
lines in the near-UV and visible spectral regions, respec- 
tively. The log gf values of 139 Fe II lines were updated ac- 
cording to the normalization given in Melendez & Barbuy 
(2002). Line by line fits to the solar spectrum were done 
in specific spectral regions, and in this process ~ 240 
lines had their log gfs selected among values from Sncdcn 
et al. (1996), Kupka et al. (1999; VALD), Biehl (1976), 
Steffen (1985), McWilliam & Rich (1994) and Bensby et 
al. (2003). Additionally, 21 hyperfine structures for the 
heavy elements Ba, La, Eu, Co, Pb were included, adopt- 
ing a solar isotopic mix. For the hyperfine structures, a 
code described and made available by McWilliam (1998) 
is employed, as described in Allen (2005). 



Accurate values for the damping constant 7 are also 
important for the computation of line profiles. This is cru- 
cial for the strong lines, and it also has a non-negligible 
effect on weaker lines (see Ryan 1998 for a discussion). 
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Fig. 2. Same as Fig. 1 for the region 610-615 nm. 



In our atomic line list, 36% of the lines had the colli- 
sional broadening obtained from Anstee & O'Mara (1995), 
Barklem & O'Mara (1997), Barklem, O'Mara & Ross 
(1998) and Barklem et al. 2000 (this series of papers is 
referred to hereafter as ABO). Those lines correspond to 
the totality of the strong neutral lines. The other lines 
were either of higher ionization stages 1 or were not in the 
energy domain range of the tabulated cross sections. Thus 
these lines were either manually fitted to the Sun or were 
assumed to have interaction constant C$ — 0.3 x 10 -31 . 

Around 5% of the lines which had Cq derived from 
ABO cross sections had to be manually fitted in order 
to reproduce the solar spectrum. This happened due to 
the dependence of the line profile both on the collisional 
broadening and on the model atmosphere employed. The 
conversion of the ABO cross sections into C§ values are 
detailed in Appendix A. 



Cross sections for ionized lines are not easily computed 
(see Barklem & O'Mara 1998 and Barklem et al. 2000). The 
cross sections for ionized transitions presented in Barklem et 
al. (2000) were included in our line list. 



The high quality of the final line list can be assessed 
by comparisons with high resolution spectra of the Sun 
(Kurucz et al. 1984) and Arcturus (Hinkle et al. 2000). 
The resolution of these spectra are of the order R ~ 10 5 , 
thus allowing a thorough verification of the accuracy of 
our wavelengths, oscillator strengths and damping con- 
stants. In Figures ^ an( i [21 comparisons of our synthetic 
spectra to those of the Sun and Arcturus are shown in 
two spectral regions. For the Sun the synthetic spectrum 
was computed adopting T e s = 5770K, log g = 4.44, mi- 
croturbulent velocity v t = 1.0 km/s and abundances from 
Grevesse & Sauval (1998). For Arcturus the parameters 
adopted were T e ff = 4275 K, log g = 1.55 and v t = 1.65 
km/s (Melendez et al. 2003). An ATLAS9 model atmo- 
sphere with these parameters was kindly computed by F. 
Castelli (private comunication) . For the abundances, we 
adopted [Fe/H] = -0.54, [C/Fe] = -0.08, [N/Fe] = +0.30, 
[O/Fe] = 0.43 and [Ni/Fe] = +0.02 from Melendez et 
al. (2003), [Mg/Fe] = +0.30 and [Ca/Fe] = +0.17 from 
McWilliam & Rich (1994). 

In both figures the Sun is very well reproduced. For the 
case of Arcturus, the synthetic spectrum looks somewhat 
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Fig. 3. Each of the seven TiO systems considered in 
this work, computed separately. The identification of each 
band is given in the figure. These spectra were computed 
with (T^, log g, [Fe/H], [a/Fe]) = (3500, 0.0, 0.0, 0.0). 

strong-lined. At such a high resolution, small errors in 
the atmospheric parameters can easily account for this 
difference, and the overall agreement is very good. 

When comparing synthetic spectrum against observed 
spectrum, one also has to be aware that observed spec- 
tra taken through the atmosphere include telluric lines, 
whereas computed spectra do not. A crucial aspect of our 
procedure is the derivation of reliable oscillator strengths 
from comparison with ground-based observations of the 
Sun and Arcturus, as explained in Barbuy et al. (2003). 
This however is not feasible in regions of the spectra which 
are heavily contaminated by telluric lines. For that reason, 
there is a gap in the line list between 1.35 and 1.50 /im - 
where the telluric absorption is severe - since we refrained 
from adopting lines with purely theoretical log gfs in that 
region. 

2.2. Molecular line opacities 

Molecular lines are an important source of opacity in the 
atmospheres of late-type stars. Molecular bands in fact 
dictate the overall shape of the spectra of late-K and M 
stars, which are the dominant contributors to the inte- 
grated light of old and intermediate-age stellar popula- 
tions longward of ~ 650 nm. Therefore, it is crucial to 
model their intensities correctly if one wants to make ac- 
curate predictions in these spectral regions. The molecules 
considered in the calculation of the present library are 
listed in Table [2 

The line list for the molecules MgH, C2 and CN-red 
were adopted from Balfour & Cartwright (1976), Phillips 
& Davis (1968) and Davis & Philips (1963), respectively. 
For the CN-blue, NH and OH-blue systems, the line 
lists by Kurucz (1993) were adopted, as implemented 
in Castilho et al. (1999). For the CH blue systems, the 
LIFBASE program of Luque & Crosley (1999) was used. 



The OH vibration-rotation line list is the one implemented 
in Melendez, Barbuy & Spite (2001). It is based on the 
laboratory work of Abrams et al. (1994), complemented 
with the theoretical line list by Goldman et al. (1998). 
The molecular gf-values of the OH lines was obtained 
from the Einstein coefficients calculated by Goldman et 
al. (1998). The CN and CO infrared line lists were im- 
plemented in Melendez & Barbuy (1999). The CO lines 
are from the theoretical work of Goorvitch (1994) and 
the infrared CN lines are from S. P. Davis, based on the 
analyses by M. L. P. Rao (S. P. Davis, private communi- 
cation). The FeH (A 4 A-X 4 A) line list of rotational lines 
were adopted from Philips et al. (1987) as described in 
Schiavon et al. (1999). The dissociation potential values 
adopted for all molecules are reported in Table 1, together 
with corresponding sources. 

The line list for the TiO 7 system was adopted from the 
theoretical work by Jorgensen (1994). The electronic os- 
cillator strength was calibrated to the intensity of the lab- 
oratory line list used in Erdelyi-Mendes & Barbuy (1989). 
For the remaining TiO electronic systems, the list by Plez 
(1998) was adopted. The Plez line list contains several mil- 
lion lines which makes the computations extremely time- 
consuming so that we decided to eliminate very weak lines 
in order to speed up our calculations. Therefore, we in- 
clude only the lines which satisfy simultaneously the fol- 
lowing set of conditions: log gf > -4.0, J" < 120, and 2/, 
v" < 9, where J" is the rotational quantum number of 
the lower level of the transition and v' , and v" are the vi- 
brational quantum numbers of the upper and lower levels 
of the transition. These values were chosen after several 
tests were performed in order to make sure that removal of 
lines not satisfying these criteria would have a negligible 
impact on our computations. 

The electronic oscillator strengths (f e i) for the TiO 
systems were obtained empirically by matching band 
strengths in spectra of observed cool stars with our syn- 
thetic spectra. For this purpose, cool giants covering the 
effective temperature range 3000 < T ff < 3800 K were 
selected from the STELIB and Indo-US libraries. In this 
temperature range, TiO bands dominate the line opacity 
in the visible region of the spectra of M giants. Because 
the model atmospheres by Castelli & Kurucz (2003) do 
not go below T c g = 3500K, we employed MARCS model 
atmospheres from Plez et al. (1992) to calibrate these f e i. 

In Figure [3] the relative contribution of each of the 
TiO systems, calculated for a cool giant, is presented. 
Some comparisons between empirically calibrated TiO- 
dominated synthetic spectra to the spectra of some ob- 
served stars are presented in Figures 0] to Moreover, the 
comparison to the spectra of an F dwarf and a K giant 
are presented in Figures [S] and [5] respectively. The wave- 
length range plotted in the later two figures were chosen 
to illustrate the spectral region where those spectral types 
dominate the flux of integrated population spectra. The 
synthetic spectra shown in these figures are the ones whose 
atmospheric parameters are the closest to those of the ob- 
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Fig. 4. Comparison between the star HD 102212 (thin line), atmospheric parameters (T e ff, log g, [Fe/H]) = (3660, 
0.0, 0.0), obtained from STELIB library and a synthetic spectrum (thick line), parameters (T e ff, log g, [Fe/H], [a/Fe]) 
= (3600, 0.0, 0.0, 0.0). This wavelength range illustrates the TiO 7 and 7' systems. 



served stars, as given in the respective library (STELIB 
or Indo-US). 

In all figures from 01 to ED the synthetic spectra were 
convolved to match the different resolutions of the ob- 
served spectra (FWHM w 1 A for Indo-US and sa 3 A for 
STELIB). 

3. Library of synthetic spectra 

Given the atomic and molecular line lists described in the 
previous section, a library of synthetic spectra was com- 
puted from 3000 A to 1.8 /Ltm, in steps of 0.02A. The model 
atmosphere grid is available for microturbulent velocities 
v t = 2 km/s, but the synthetic spectra were calculated 
with values more representatives of observed stars as given 
below: 

— Vt = 1.0 km s" 1 for log g > 3.0; 

— Vt = 1.8 km s _1 for 1.5 < log g < 2.5, and; 

— Vt = 2.5 km s _1 for log g < 1.0. 

The output of the synthesis code is normalized flux. 
In order to produce spectra in an absolute flux scale, the 
normalized high-resolution spectra were multiplied by the 
true continuum given directly by the ATLAS9 model at- 
mosphere. The library is presented in both normalized 
and absolute flux formats, and each spectrum is stored 
in a FITS file consisting of two apertures (multispec IRAF 



format), where aperture number 1 is the normalized spec- 
trum and aperture number 2 is an absolute flux spectrum. 
The spectral library is available upon request. 

3.1. Coverage in Stellar Parameter Space 

The main grid covers the following parameters: 

— Effective temperatures: 3500 < T e g < 7000 K in steps 
of 250K 

— Surface gravities: 0.0 < log g < 5.0 in steps of 0.5 

— Metallicities: [Fc/H] = -2.5, -2.0, -1.5, -1.0, -0.5, 0.0, 
0.2 and 0.5 

— Chemical compositions: [a/Fe] = 0.0 and 0.4, where 
the a-elements considered are O, Ne, Mg, Si, S, Ca 
and Ti. 

It is not well established yet whether titanium is as- 
sociated to the synthesis of the a-elements or to the iron- 
peak elements (Timmes et al. 1995). Nevertheless, it has 
been shown that its abundance follows the same pattern 
as that of a-elements for several types of populations (e.g. 
Mc William & Rich 1994; Pompeiaet al. 2003), and thus we 
decided to include this element in the a-cnhanced group. 

Concerning calcium, observations of halo and disk 
Galactic stars show that [Ca/Fe] follows closely the be- 
haviour of other a-elements as a function of [Fe/H] 
(Mc William & Rich 1994; Edvardsson et al. 1993). In the 
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Fig. 5. Comparison between the star HD 39801 (thin line), atmospheric parameters (T ff, log g, [Fe/H]) = (3540, 1.0, 
0.0), obtained from Indo-US library and a synthetic spectrum (thick line), parameters (T e g, log g, [Fe/H], [a/Fe] 
(3600, 1.0, 0.0, 0.0). This wavelength range illustrates the TiO a system. 



bulge of our galaxy, however, some high-resolution anal- 
yses show that [Ca/Fe] seems to be close to solar (e.g. 
Zoccali et al. 2004). From lower resolution studies, there is 
evidence as well that [Ca/Fe] is not enhanced in elliptical 
galaxies (e.g. Thomas et al. 2003b), although Prochaska, 
Rose & Schiavon (2005) showed that the calcium abun- 
dance derived from the Ca4227 index is dependent on the 
index-measuring strategy. 

Considering the application of this grid to the study of 
old and metal-rich stellar populations, a subsample of the 
a-enhanced main grid was also computed with [Ca/Fe] = 
0.0, for [Fe/H] = -0.5, 0.0, 0.2 and 0.5. The sets of T off and 
log g for these latter computations were chosen amongst 
those that are relevant for old stellar populations at these 
metallicities, for ages between 8 and 13 Gyr. Bertelli et al. 
(1994) isochrones were used to guide our selection. Some 
of these isochrones are presented in Figure UHl as well as 
the (T e ff, log g) pairs computed for this extension of the 
library. 

High resolution synthetic spectra computed for a cool 
giant with solar-scaled and a-enhanced abundance pat- 
terns are compared in II II and 1121 for the metal-poor and 
metal-rich cases, respectively. As expected, well known 
indicators of abundances of a-elements, such as the Mg 
triplet lines at 517 nm, the MgH band at 518 nm and the 
Ca triplet in the far red, are stronger in the a-enhanccd 
spectra. We note however other lines which are very sensi- 
tive to a-enhancement and yet are virtually unexplored in 
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Fig. 10. Bertelli et al. (1994) isochrones for three metal- 
licities ([Fe/H] = -0.4, 0.0, +0.4; solid, dotted and dashed 
lines) and three ages (8, 10.5 and 13 Gyr; from left 
to right). The circles indicate the extension of the a- 
enhanced main grid, computed with [a/Fe] = 0.4, [Ca/Fe] 
= 0.0 and [Fe/H] = (-0.5, 0.0, 0.2 and 0.5). 
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Fig. 6. Comparison between the star HD 126327 (thin line), atmospheric parameters (T e ff, log g, [Fe/H]) = (3000, 0.0, 
-0.58), obtained from Indo-US library and a synthetic spectrum (thick line), parameters (T e ff, log g, [Fe/H], [a/Fe]) 
= (3000, 0.0, 0.0, 0.0). This wavelength range illustrates the TiO a system. 



stellar populations work (e.g. the Ca H & K lines and the 
lines due to the CN molecule, at A > 1093 nm). We also 
note that some lines are visibly weaker in the a-enhanced 
case (e.g. the region around 385 nm in Figure EJ- This 
cannot be attributed to individual abundance variations 
of the non-a elements, because those abundances are the 
same as in the solar scaled case. The decrease of the inten- 
sity of some lines can be explained because the a-elements, 
specially Mg, are electron donors and thus contribute in 
an important way to the continuum opacity. Because the 
spectra shown in Figures 1111 and 1121 are normalized, the 
net effect is that some lines in the a-enhanced case look 
weaker due to enhanced continuum opacity, even when the 
individual abundances corresponding to those lines are un- 
changed. 



3.2. The relative contribution of different chemical 
species to the total line opacity 

A great advantange of synthetic spectra over observed 
spectra, is that with spectral synthesis it is possible to 
disentangle the relative contribution of different chemical 
species to the line spectrum, which allows the study of 
line indices in a more careful way. For example, it is pos- 
sible to evaluate the presence of blends or the range of 
applicability of a given line indice. 



As an example, the well-known Lick indices Mgi , Mg2 
and Mgh measure the intensity of the Mgi triplet lines 
(5167.327, 5172.698 and 5183.619 A) together with MgH, 
C2 and TiO bands, besides other atomic lines. These in- 
dices measure essentially Mgi and MgH in K giants and 
in the integrated light of galaxies. However, it should be 
noted that in stars cooler than 3750 K, bands from the 
a system of TiO become stronger than MgH, as can be 
seen in Figure El It is clear that for very metal-rich stars 
at the tip of the red giant branch and galaxies containing 
such stars these indices are substantially affected by TiO 
bands, which makes their interpretation somewhat more 
complicated (see also Schiavon 2005). 

The separate contribution of atoms and molecules is 
shown in Figure El for the parameters (T e g , log g) = 
(3500, 0.5), (4750, 2.5) and (6000, 4.0) for [Fe/H] = -2.0 
and 0.0 ([a/Fe] = 0.0). In order to illustrate the contribu- 
tion of atoms and each diatomic molecule, we computed 
each molecule separately. The results are shown in Figures 
[TBI - El where the impact of each molecule on the line 
spectrum can be evaluated. 

4. Summary 

A library of synthetic stellar spectra is presented. This 
library was computed with a sampling 0.02A and for the 
parameters: effective temperatures 3500 < T e g < 7000 
K, surface gravities 0.0 < log g < 5.0, metallicities -2.5 
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Fig. 7. Comparison between the star HD 217906 (thin line), atmospheric parameters (T e ff, log g, [Fe/H]) = (3600, 1.2, 
-0.11), obtained from Indo-US library and a synthetic spectrum (thick line), parameters (T e ff, log g, [Fe/H], [a/Fe]) 
= (3600, 1.0, 0.0, 0.0). This wavelength range illustrates basically the TiO (3 system. 



< [Fe/H] < +0.5, and a-elements over Fe ratio [a/Fe] = 
0.0 and 0.4. For the metallicites [Fe/H] > -0.5, two en- 
enhanced grids are available, with and without Ca in the 
enhanced group. 

For the computation of this grid, an atomic line list was 
implemented based on the previous lists of Barbuy et al. 
(2003) and Melendez & Barbuy (1999). The log gf values 
were updated to the most recent NIST values, and damp- 
ing constants were calculated for the totality of strong 
lines following ABO papers. 

The molecular line opacities include the molecules 
MgH, C 2 , CN, CH, CO, OH, NH, FeH and TiO. The TiO 
systems, which dominate the opacity of cool stars in the 
visible wavelength range, were calibrated by comparison 
with observed stars in the empirical libraries STELIB and 
Indo-US. 

This high-resolution stellar synthetic library covers si- 
multaneosly a wide wavelength range, and a-enhanccd 
chemical compositions. The a-enhancement in particu- 
lar is indispensable to the reproduction of the integrated 
properties of old stellar populations, for all metallicities 
regimes. This library is a valuable and powerfull tool for 
the study of stellar populations, and can be used to model 
old metal-rich stellar populations with unprecedented ac- 
curacy. 

Note added after acceptance: We would like to mention 
that by the time of the acceptance of our paper, two recent 
libraries including variable a-enhancement came to our 



knowledge, Munari et al. (2005) and Brott 
(2005). 
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Appendix A: Calculation of collisional broadening 

The main source of collisional broadening for metallic lines 
in stars of spectral types between F and M is collisions 
with neutral hydrogen. The spectral synthesis code eval- 
uates 7 for each atmospheric layer given the interaction 
constant Cq: 



(A.l) 



where v is the relative velocity between the collid- 
ing particles and Nh is the density of Hydrogen atoms. 
The Cq parameter (or any other parameter related to the 
damping width) is accurately known only for very few 
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Table 1. Molecules considered in the computations of the synthetic spectra. References: (1) Huber & Herzberg (1979); 
(2) Pradhan et al. (1994); (3) Schulz & Armentrout (1991); 



Molecule 


System 


Number 


Wavelength 


Number of 


Dissociation potencial (ref) 






of lines 


coverage (nm) 


vibrational bands 


Do (eV) 


MgH 


A 2 n-X 2 E 


1945 


404 - 609 


13 


1.34 (1) 


C2 Swan 


A 3 n-x 3 n 


11254 


429 - 676 


35 


6.21 (1) 


CN blue 


B 2 E-X 2 E 


92851 


300 - 600 


197 


7.72 (2) 


CN red 


A 2 n-X 2 E 


23828 


404 - 2714 


74 




CH AX 


A 2 A-X 2 n 


10137 


321 - 786 


20 


3.46 (1) 


CH BX 


B 2 A-X 2 n 


2016 


361 - 682 


10 


J) 


CH CX 


c 2 A-x 2 n 


2829 


270 - 424 


12 


)} 


CO nir 


X X E+ 


7088 


1578 - 5701 


63 


11.09 (1) 


NH blue 


A 3 n-X 3 E 


8599 


300 - 600 


55 


3.47 (1) 


OH blue 


A 2 E-X 2 n 


6018 


300 - 540 


46 


4.39 (1) 


OH nir 


x 2 n 


2028 


746 - 2594 


43 




FeH 


A 4 A-X 4 A 


2705 


778 - 1634 


9 


1.63 (3) 


TiO 7 


A 3 $-X 3 A 


26007 


622 - 878 


23 


6.87 (1) 


TiO 7' 


B 3 n-X 3 A 


219367 


501 - 915 


81 




TiO a 


C 3 A-X 3 A 


360725 


380 - 861 


79 




TiO p 


c 1 $-a 1 n 


91804 


431 - 804 


63 




TiO S 


b 1 n-a 1 A 


189019 


622 - 1480 


66 




TiO e 


E 3 n-X 3 A 


253755 


641 - 1341 


61 




TiO <f> 


^n-tPs 


105082 


665 - 1780 


65 






I I I I I I I I I I I 

400 450 500 

Wavelength (nm) 



Fig. 8. Comparison between the star HD 126141 (thin line), spectral type F5V and atmospheric parameters (TW, log 
g, [Fe/H]) = (6670, 4.3, -0.02), obtained from STELIB library and a synthetic spectrum (thick line), parameters (T e ff, 
log g, [Fe/H], [a/Fe]) = (6750, 4.5, 0.0, 0.0). 
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Fig. 9. Comparison between the star HD 148513 (thin line), spectral type K4III and atmospheric parameters (T ff, 
log g, [Fe/H]) = (3979, 1.03, -0.14), obtained from STELIB library and a synthetic spectrum (thick line), parameters 
(T cff , log g, [Fe/H], [a/Fe]) = (4000, 1.0, 0.0, 0.0). 



lines. The Unsold (1955, see also Gray 1976) approxima- 
tion allows one to straightforwardly estimate the damping 
constants, but it has long been recognized that it under- 
estimates the collisional broadening, and thus empirical 
enhancement factors are applied to Cq in order to repro- 
duce the observed line profiles. 

The works by O'Mara (1976) and Anstee & O'Mara 
(1991) developed an alternative theory for the computa- 
tion of cross sections a of atomic transitions, thus provid- 
ing a more accurate way to evaluate 7. Anstee & O'Mara 
(1995), Barklem & O'Mara (1997) and Barklem, O'Mara 
& Ross (1998) published cross sections for a wide range 
of s-p and p-s, p-d and d-p, d-f and f-d transitions, re- 
spectively. Herefrom, we refer to these cross-sections as 
ABO. 

In order to use the more precise ABO theory in our cal- 
culations, we had to relate the ABO cross section a to the 
classical interaction constant Cq which is the input value 
to our spectrum synthesis code. The interaction constant 
can be expressed by: 



where i??. and Rf a are the mean-square radii of the 
upper and lower states of the transition in atomic units. 
The cross section from the van der Waals theory is given 
by: 

a vdW = 63.65(i?L - Rj fl b (A.3) 

From ea uat ions I A . 21 and I A . 31 and replacing the classi- 
cal cross-section a v dw by the cross-section given by ABO 
theory (TabOi we obtain the relation: 

C 6 = 6.46 x 10" 34 (cr ABO /63.65) 5/2 



C 6 = 6.46 x lQ- 34 (Rl i - R 



loJ 



(A.2) 



(A.4) 

Or in terms of 76: 

l6 /N H = 1.415 x 10- U v 3 / 5 a abo (A.5) 

The expression in Ea. I A. II leads to a 7 oc T 2 / 5 relation- 
ship. However, according to the theory developed by ABO, 
which is based on the formalism by Brueckner (1971), 7 
has a dependence on the temperature that varies between 
T ' 3 and T 0A . The ABO papers tabulate the velocity ex- 
ponents a that take into account this non-universal depen- 
dence. Due to the different formalisms used by ABO the- 
ory and by the spectral synthesis code, the dependences of 
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Fig. 11. Solar scaled (thin line) and a-enhanced (thick line) synthetic spectrum computed with T e ff = 4500K, log g 
= 3.0 and [Fe/H] = -2.0, for different wavelength regions. The spectra were convolved to FWHM ~ 0.5 A. 



76 on temperature are somewhat different. Nevertheless, 
the error introduced in the computed line profile is negligi- 
ble (< 5%) given other sources of errors like uncertainties 
in the atmospheric parameters, model atmospheres and 
NLTE effects. 

The a abo values were obtained primarily from 
Barklem et al. (2000). Cross sections for lines not included 
in Barklem et al. (2000) were calculated through the code 
presented in Barklem, O'Mara & Ross (1998). For this 
purpose, the energy levels and orbital angular momentum 
quantum numbers for each transition were obtained from 
the line list of Kurucz (1993), available at the web address 
http://kurucz.harvard.edu/linelists The ionization limits 
of the transitions were approximated by the first ioniza- 
tion potential of the atom. Although this approximation 



is not strictly accurate, as pointed out in Barklem et al. 
(2000), the detailed inclusion of the ionization limits of 
each transition is beyond the scope of the present work. 

A further step in critically evaluating the damping 
widths of the strongest lines is to take into account the 
model atmosphere adopted. This effect is illustrated for 
the Ca line 6 162 A in Figure lA~Tl Cross sections for three 
lines of Ca I (A = 6102, 6122, 6162) have been accu- 
rately calculated by Spielfiedel et al (1991). The ratios 
1 abo h 'Spiel fiedei for these lines are ( 0.8, 1 .05, 0.96), which 
indicates good agreement. In Figure lA~T1 is shown the so- 
lar spectrum (solid line), a synthetic spectrum calculated 
using the solar model atmosphere by Holweger & Mullcr 
(1974, HM74, filled circles), and another synthetic spec- 
trum calculated with ATLAS9 models (open circles) . The 
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Fig. 12. Solar scaled (thin line) and a-enhanced (thick line) synthetic spectrum computed with T e ft = 4500K, log g 
= 3.0 and [Fe/H] = 0.0, for different wavelength regions. The spectra were convolved to FWHM ~ 0.5 A. 



Ce employed is the one directly obtained from equation 
I A. 41 using the same parameters as given in Anstee & 
O'Mara (1995), Table 3. 

It can be seen that the line profile computed with the 
ATLAS9 model is stronger than when the HM74 model is 
used. The use of a MARCS model atmosphere (Plez et al. 
1992) results in a line profile that is very similar to that 
obtained when the ATLAS9 model is employed. 

The computed line profile depends on both the damp- 
ing constant and the model atmosphere. Hence, even em- 
ploying perfect damping parameters does not guarantee 
that correct profiles are computed, since the synthetic pro- 
file also depends on the adopted model atmosphere, as al- 
ready pointed out by Cayrel et al. (1996) and Barbuy et 
al. (2003), and as is clearly shown in Fig. Al. 



As a result, 5% of the lines for which the broaden- 
ing values were calculated through ABO theory have too 
strong wings when compared with the observed spectrum 
of the Sun, when calculated with an ATLAS9 model. 
Therefore, the C§ values for these lines were fitted manu- 
ally to the solar spectrum. It was found that the average 
7(ABO) /7(best fit) for those lines is « 1.4. Cayrel & Van't 
Veer (private communication) performed similar tests for 
the Mg I triplet lines, and found that 7(ABO)/7(best fit) 
sa 1.5, which is in remarkable agreement with our own 
results. 

The wings of the weaker lines are not critically affected 
by the use of different model atmospheres (as can be seen 
for the two weaker lines in Figure lA~T|) . This implies that 
for w 95% of the lines, the broadening predicted by ABO 
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Fig. 13. Individual contribution of the MgH and TiO a molecules, for three T e ff and two log g (indicated in the 
figures). The top plot in each panel represents the MgH spectrum, and the bottom plot presents the TiO a spectrum. 



theory fits well the solar spectrum when ATLAS9 models 
are employed. 
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have non-negligible effect in the total spectrum. 



Edited by A. McWilliam and M. Rauch. Pasadena: 

Carnegie Observatories 
Melendez, J., Barbuy, B., 1999, ApJS, 124, 527 
Melendez, J., Barbuy, B., 2002, ApJ, 575, 474 
Melendez, J., Barbuy, B. & Spite, F. 2001, ApJ, 556, 858 
Melendez, J., Barbuy, B., Bica, E., Zoccali, M., Ortolani, S., 

Renzini, A., Hill, V., 2003, A&A, 411, 417 
Mendes de Oliveira, C, Coelho, P., Gonzalez, J. J., Barbuy, B., 

2005, AJ, 130, 55 
Mu nari, U., Castelli, F., Zwitter, T., 2005, A&A, accepted 

( astro-ph/0502047 ) 
Murphy, T., Meiksin, A., 2004, MNRAS, 351, 1430 
O'Mara, B., 1976, MNRAS, 177, 551 

Panter, B., Heavens, A., Jimenez, R., 2004, MNRAS, 355, 764 
Pickles, A.J. 1998, PASP, 110, 863 
Plez, B., 1998, A&A, 337, 495 



Plez, B., Brett, J.M., Nordlund, 1992, A&A, 256, 551 
Phillips, J. G., Davis, S. P., 1968, Berkeley Analyses of 

Molecular Spectra, Berkeley: University of California Press 
Phillips, J.G., Davis, S.P., Lindgren, B., Balfour, W.J. 1987, 

ApJS, 65, 721 

Pompeia, L., Barbuy, B., Grenon, M., 2003, ApJ, 592, 1173 
Praderie, F., 1967, Annales d'Astrophysique, Vol. 30, p. 31 
Pradhan, A.D., Partridge, H, Bauschlicher, C.W.J. 1994, J. 

Chem. Phys., 101, 3857 
Prochaska, L.C., Rose, J. A. & Schiavon, R.P., 2005, AJ, sub- 
mitted 

Proctor, R. N., Forbes, D. A., Beasley, M. A., 2004, MNRAS 
355, 1327 

Reader, J., Wiese, W. L., Martin, W. C, Musgrove, A., Fuhr, J. 
R., 2002, NASA Laboratory Astrophysics Workshop, held 
May 1-3 2002 at NASA Ames Research Center, Moffett 



P. Coelho et al.: A library of synthetic spectra with a-enhancement 



19 




616.1 616.15 616.2 616.25 616.3 

Wavelength (nm) 

Fig. A.l. Comparison between the solar spectrum (solid 
line, no symbols) and synthetic spectrum calculated with 
ATLAS9 (solid line, open circles) and HM74 (solid line, 
filled symbols) model atmospheres, around the Ca 6162 A 
line. The same parameters were adopted for the atomic 
transition considered in both calculations, and thus the 
difference in the line profiles is due only to the different 
model atmospheres employed. 
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